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The knowledge of a series expansion of the equation of state provides a deep insight into the 
physical nature of a quantum system. Starting from a generic "perturbative" equation of state of 
a homogeneous ultracold gas we make predictions for the properties of the gas in the presence of 
harmonic confinement. The local density approximation is used to obtain the chemical potential, 
total and release energies, Thomas- Fermi size and density profile of a trapped system in three-, two-, 
and one- dimensional geometries. The frequencies of the lowest breathing modes are calculated using 
scaling and sum-rule approaches and could be used in an experiment as a high precision tool for 
obtaining the expansion terms of the equation of state. The derived formalism is applied to dilute 
Bose and Fermi gases in different dimensions and to integrable one-dimensional models. Physical 
meaning of expansion terms in a number of systems is discussed. 



I. INTRODUCTION 



Recent progress in experimental techniques has made it possible to investigate the properties of gases of bosonic 
atoms in traps at temperatures lower than that required to achieve Bose-Einstein condensation (BEC), as well as the 
properties of fermionic atoms at temperatures lower than the Fermi temperature. When the density of a homogeneous 
gas is small enough perturbation theory can be applied and the equation of state of repulsive bosons and two- 
component fermions can be expressed in terms of the expansion parameter naf, where n is the density and a s is the 
s-wave scattering length. The problem of a cold trapped gas is more difficult as another parameter, the harmonic 
oscillator length a/j , enters into the problem. Although many properties of a homogeneous system (equation of state, 
etc.) can be predicted quite easily, no exact solutions in the presence of an external potential are generally known. 
If the number of particles is sufficiently large and the energy per particle is much larger than the level spacing of the 
harmonic oscillator, one can apply the local density approximation (LDA). This method is commonly and successfully 
used to relate properties of an inhomogeneous system to the ones of the corresponding homogeneous system. In this 
way the predictions of known homogeneous theory can be tested experimentally in confined systems. 

Alternatively, the problem can be seen from the other side. It might happen that it is possible to observe and 
investigate properties of a trapped system in an experiment, while an analytical description of the same physical 
system is not known even for a homogeneous case. Then the experimental study of a trapped system might shed 
light on the properties of a homogeneous system and allow us to judge the validity of the analytical approaches. For 
example, this is the case with the two-component Fermi gas in the BEC-BCS crossover. The unitarity regime was 
already observed in experiments (see, for example, while there is a lack of a well-posed and complete theory. 

Another intriguing question is presence or absence of the Lee- Yang term, corresponding to a beyond mean field term 
of composite bosons, in the equation of state on the BEC side where formation of molecules was observed H,B,|^. 
This problem is still an open question and there is no conclusive theoretical answer. A high precision experiment 
might resolve the question. 

The aim of the current work is to investigate how the expansion of the equation of state of a homogeneous system 
shows up in properties of a trapped system. In order to do so we assume that the equation of state of a homogeneous 
gas can be written as as series expansion in some limit. The local density approximation is used to obtain measurable 
energetic and spatial properties of a trapped system as a function of the parameters of the expansion and LDA 
parameter. Special attention is paid to the prediction of the frequencies of collective oscillations, as a sensitive 
experimental tool for investigation of the equation of state. Scaling and sum-rule approaches are used for deriving 
the frequencies of small collective oscillations. The derived formalism is applied to dilute Bose and Fermi gases in 
different dimensions and to integrable one-dimensional models. We discuss the physical meaning of expansion terms 
in a number of systems. 

The paper is organized as follows. In Section [n] we explain the local density approximation in an arbitrary number 
of dimensions and introduce dimensionless units. In Sec. IIIII we introduce the "perturbative" equation of state and 
solve the local density approximation problem. Explicit expressions for the chemical potential, energy, release energy, 
size and density profile of a trapped gas are obtained. In Sec. IIVI we use scaling and sum-rule approaches for the 
prediction of the frequencies of collective oscillations. Cases of spherical, elongated, "cigar"- and "pancake"- shaped 
traps are considered. Formulas are derived for three-, two- and one- dimensional traps. The LDA is applied to 
relate the "perturbative" equation of state to frequencies of the breathing mode. In Sec. we apply the developed 
formalism to a number of physical systems where perturbation expansion of the equation of state is known. These 
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include dilute bosonic gases (in 3D, quasi 2D, and quasi ID geometries), two-component Fermi gases (in 3D and in 
ID), as well as a number of integrable one-dimensional models (Lieb-Liniger, Calogero-Sutherland). In Sec. I VII we 
draw our conclusions. 



II. LOCAL DENSITY APPROXIMATION 



We consider a gas of N particles (bosons or fermions) at a density n = N/V. We suppose that the density 
dependence of the chemical potential fihom of the homogeneous system is known from either perturbation theory or 
from a fit to experimental data. Within the local density approximation one assumes (see, for example, |(j) that the 
chemical potential /i of a trapped system is given by the sum of the local chemical potential, which is taken to be 
equal to the chemical potential of the uniform system fXhom at the corresponding density, and the external field V ex t- 



A* = Mfcom(«(r)) + V ext (r) 



(1) 



In the following we will always consider a harmonic external confinement for two basic reasons. First, this is the 
most widely used trapping potential in cold gases. Second, a quadratic potential can be used to describe the external 
potential close to its minimum. In order to consider the most general case we will assume that the trapping frequencies 
are different, so that the confinement is anisotropic V ex t(r) — ^mui x x 2 + ^mLu y y 2 + ^mcu z z 2 . 

The value of the chemical potential /i of a trapped gas is fixed by the normalization condition 



N = I n(v)dr = 
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Here the density profile is obtained by inverting the density dependence of the local chemical potential n = n^ oia . 

The spatial extensions of the cloud, R x , R y , R z , depend on the value of the chemical potential and frequencies of 
the harmonic confinement: 



/i = ^mujlR 2 x 
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(3) 



By expressing the distances in the trap in units of the size of the cloud: f — (x / R Xl y / R v , z / R z ) it follows imme- 
diately that the frequencies of the trap enter in Eq. J3J) only through the combination uj) lo — (uJxUjyUJz) 1 /^ and the 
normalization condition can be conveniently written in terms of the averaged oscillator length a/ lo = yjh/mujho- Now 
the integral is to be taken inside a sphere of radius 1 and is symmetric with respect to f. 

The normalization condition (J2J) can be conveniently rewritten in dimcnsionlcss units: 



k D D ' 2 = li D ' 2 



J /C - r 2 )] dr, n(r) = for |r| > 1 



(4) 



Here the chemical potential p, is expressed in units of ^N 1 ' D huhoi which is a natural unit of energy in the trap, 
the density in a homogeneous system is measured in units of a~ D , where a is a length scale convenient for the 
homogeneous system (for example, it can be equal to the s-wave scattering length a s ), chemical potential (i.e. the 
argument of the inverse function fJ^ m ) is measured in units of Ti 2 /ma 2 , D — 1, 2, 3 is the dimensionality of the system, 
and the LDA characteristic parameter Ad is defined as 



A 



D 



(•5) 



From the Eq. Q it becomes clear that within the LDA there is a scaling in terms of the characteristic parameter A^. 
In other words systems having a different number of particles and different frequencies of the harmonic confinement 
will have the same density profile and other LDA properties (once expressed in appropriate units as discussed above) 
if they have equal values of the parameter A^. 



III. SOLUTION FOR "PERTURB ATIVE" EQUATION OF STATE 

In this Section we develop a general theory for a "perturbative" equation of state in a three-, two-, and one- 
dimensional systems. We start by writing down the equation of state /Jhom(")- Assuming that a perturbation 
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approach can be applied to a homogeneous system, we suppose that /ihom(^) is known as a series expansion, either at 
small or at high density, in terms of the expansion parameter (it can be gas parameter naP in the weakly interacting 
regime and l/(na D ) in the regime of strong interactions): 



Mhom = C x (na D p(l + C 2 {na D y* + ...)- 



(6) 



In this formula the second term is considered to be small compared to the first term (|C2(na 3 ) 72 1 <C 1) and will be 
treated perturbatively in the following. We will limit the series to the second term, as it is responsible for the main 
contribution to the dependence of the frequency of the breathing mode on the density, and we omit higher terms in 
the expansion (BJ. 

We obtain the value of the chemical potential of a trapped system by resolving the normalization condition |0J. 
The result can be conveniently expressed in terms of the characteristic parameter A/j (|5l). 
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where constants C M i, C M 2 are independent of the parameters of the trap and, instead, are directly related to the 
equation of state of a homogeneous system. In order to find them we evaluate the integral (@J in three-, two-, and 
one- dimensional geometries. The explicit expressions for the coefficients C M i and C M 2 are given in Tabic [I] The 
leading term in (0) depends only on two parameters, C\ and 71, of the homogeneous equation of state, while in order 
to obtain the next term one has to know four parameters: C\, 71, C2, 72- This type of behaviour is quite general and 
applies to majority of the quantities we estimate in a trapped gas, with the exception of the frequencies of collective 
oscillations, where the leading term depends only on 71 (see below). 
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TABLE I: Coefficients C M i and C M 2 of the expansion of the chemical potential /1 of a trapped system Q obtained from 
evaluation of the normalization integral 0. In a 3D system the integral is to be taken inside a sphere and can be evaluated 
explicitly if condition (1 +72V71 + 1 > is satisfied (see, for example, 7]). In a two-dimensional system, D — 2, the integral is 
to be calculated inside a circle of unit radius. There is a substitution f /(l — f 2 )27rf df = J f(y)^ dy which together with a 
particular form of the chemical potential ||SJ leads to an elementary integral. In a one-dimensional system, D — 1, the spatial 
integration over the coordinate z goes from —7? to R. In contrast to the two-dimensional case no simplification is found and 
the final expression for the chemical potential fi 10 will contain the gamma functions as in fully three-dimensional case. The 
integration was carried out under conditions 71 + 72 + 1 > 0. 



The energy of the trapped system can be obtained by integrating the chemical potential E = L (i(N')dN' . From 
Eqs. J2J and Q we derive an explicit expression 
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The density profile of a trapped gas can be accessed experimentally by using an absorption imaging technique, 
when the shadow of the cloud of the cloud is imaged on a CCD camera. The size of the cloud can be extracted from 
such images. In our approach the size of the cloud along axis a (a = x, y, z) is given by R a = RoJho/^a, where the 
expression for the Thomas-Fermi radius R is 



R = N^C 1 J 1 2 kl D -< 1+2 
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The applicability of the local density approximation is that the energy of the gas © is large compared to the 
interlevel spacing E 3> ^toho, so that the discretization of the levels in the trap can be ignored. The same condition it 
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terms of length means that the size of the gas © is large compared to the oscillator length R 3> aha- These conditions 
are usually satisfied for a large number of atoms. In a one-dimensional system the LDA works well even for a small 
number of atoms 0] , as effects of interactions are enhanced in low-dimensional systems compared to three-dimensional 
systems. 

The density profile n(r) has a typical shape for the local density approximation: 
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where the chemical potential is given by (JJJ and the size of the cloud is as in @. For example, for an expansion in 
the mean-field regime 71 = 1, the leading term has the shape of an inverted parabola. In the case of a ID weakly 
interacting Fermi gas or Tonks- Girardeau gas, 71 = 2, and it takes shape of a semicircle. 

Another important spatial quantity is the mean square displacement (r 2 ). It is directly related to the potential 
energy of the harmonic confinement. The integration of r 2 over the density profile (|10f> can be performed and the 
result of such an integration is conveniently written in units of the size of the cloud: 
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Another relevant quantity that can be accessed in experiments is the release energy defined as the difference of the 
total energy E of the trapped gas and the potential energy Eho = Y muJ ho( r2 ) 01 * ne confinement: E re i — E — E% . 
The potential energy of the confinement can be eliminated by switching off the trap and the energy of the expanding 
cloud is given by E re i. Combining together Eqs. iJSJ and 111|) and after some algebra we obtain the following expression 
for the dimensionless release energy E re \ = E re \ / '\N~b Tiujho'- 
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IV. FREQUENCIES OF COLLECTIVE OSCILLATIONS 



A very important aspect of the LDA approach is that it can be used for the prediction of the frequencies of 
the collective oscillations. If oscillations are generated by displacing the center of the trap, the frequency of the 
resulting excitation is determined only by the harmonic confinement and is independent of the interactions between 
atoms. Thus, this type of excitation can not be used to probe the equation of state, but helps to measure the trap's 
frequencies. Instead, if the collective mode is excited by changing slightly the frequency of the confinement, the 
resulting "breathing" oscillation depends on the speed of sound in the gas c = y/n/m d/i/dn and probes directly the 
equation of state. 

If the trap is spherical, then the mean square displacement (r 2 ) can be used to estimate the frequencies of the 
collective modes. An important point is that the breathing mode is naturally excited by a closely related operator 
'Ylii r 'i- I n this case the sum rules approach can be applied leading to a compact expression for the frequency of the 
breathing mode 
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This formula was obtained in [9| while analyzing collective frequencies of a one-dimensional Lieb-Liniger system. 
The applicability of formula (jl3|l is not restricted to a one-dimensional system and remains valid in two- and three- 
dimensional systems. 

Calculation of the derivative of the mean square displacement Ijlll) with respect to the square of the trap frequency 
leads to the following expression for the frequency Q of the breathing mode in a spherical trap: 
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The constant term in the expression for the frequencies Ijl4(l depends on the power of the leading contribution to 
the homogeneous chemical potential 10 and on the dimensionality of the system. The dependence on parameter A^ 
enters only in the subleading term which we are able to obtain with this perturbative calculation. 
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An alternative method is needed in order to study collective modes in anisotropic traps where the formula l|13[) is not 
applicable. As proved in ^(jj the scaling approach allows us to calculate the frequencies of the collective oscillations 
with extremely high precision in anisotropic harmonic confinement. The method is based on using a scaling ansatz 
for the time-evolution of the density profile n(x,y,z,t) oc n(ai(t)x, a2(t)y, az(t)z) and solving the hydrodynamic 
equations under the assumption that the amplitude of oscillations is small. The oscillation frequencies w in a 3D 
elongated anisotropic trap with trap frequencies (wj_, u>±, uj z ) are given by 
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2EluI + (2 + E)uj 2 z ± ^(25^ + (2 + 5)cj 2 ) 2 + 16(2 - 3E)u 2 ± uj 2 



(15) 



The natural parameter 5 is related to the average of the compressibility mc 2 = n$/ihom/<9n and the potential energy 
of the harmonic oscillator: S = D (tog 2 ) / (Eho) [lCL Hlj. We evaluate the parameter 5 explicitly: 

— 2(1 + 71)+ rtlMufT^j teJ A - (16) 

From Eq. I|15|l it follows that the solution in a a spherical trap lu±_ = lo z = ojho can be simplified. There are two 
frequencies: ^3£>/w 2 lo = 2 (center of the mass mode) and ^In/^ho = 3/2 5 — 1 ("breathing" mode). Direct calculation 
shows that the sum rules prediction for the spherical trap, Eq. (|14H with D = 3, coincides with the predictions of 
Eq. (|15l) . In an elongated trap, lu z + 1 it is possible to excite axial f2 z and radial f2j_ modes separately. In a "cigar" - 
shaped trap uj z <C the frequencies of these modes are ^ 2 /w 2 = 3 — 2/5 and fl]_/uj^_ = 5 correspondingly. Another 
possible configuration is a "pancake" trap, where lo z 3> u>j_. In this situation the frequencies are fl 2 /uj 2 = 1 + 5/2 
and 0^/^ = 6- 16/(2 + 5). 

We can generalize this scaling method to systems with low dimensionality. We find that in a two-dimensional 
system the frequencies of small oscillations are given by 



SI 



2 

2D 



1 



(2 + S)(u,i + ul) ± ^(2 + S) 2 K + c 2 ) 2 - 325u; 2 W 2 



(17) 



In a spherical trap the frequencies are equal to: O^-oA'ho = 2 (center of mass mode) and Ql D /u 2 = 5 ("breathing" 
mode). Again, the frequency of one type of excitation is defined only by the the trap frequency, while the frequency of 
the other mode can be used to study the equation of state of a homogeneous system. The result obtained for a spherical 
trap is the same as from Eq. (|14(l and D = 2. In an elongated trap (lu z <C uj±) one obtains the following expressions 
for the frequencies of collective excitations in terms of the parameter 5: VL 2 ] _/uj 2 l = 1 + 5/2, Vt} z jio\ — 4S/ (2 + 5). 

There is no way of having an anisotropic trap in a one-dimensional system, so the result of the scaling approach 

n 2 1D /iu 2 = 1 + 5/2 (18) 

coincides with the prediction of a spherical trap Eq. <|14(l and D = 1. 

It worth noticing that for a given equation of state © the dependence on the number of particles N, a/ lo , a (i.e. 
the interaction strength) enters only in the second term. This means that by using only the leading term in © 
one obtains no dependence on A^jj in ll(j|) and it is necessary to consider the next term to describe changes in the 
frequency. This justifies the necessity of our perturbative approach ©. 



V. DISCUSSION 



As the frequencies of the collective oscillations can be measured in experiments with high precision, formulae l|14I 
II fill provide a very useful tool for the investigation of the equation of state. As an example, the description of the 
two-component Fermi gas close to a Feshbach resonance (unitary regime) is a very non-trivial problem as there is no 
obvious small parameter that can be used to construct an analytically solvable expansion theory (see Fig. This 
regime has been studied in recent experiments [T3 . Il3l Il4| . The s-wave scattering length in the unitarity regime 
becomes larger than any other length-scale in the system \a s \ — > oo (A 3 £i — > oo) and, essentially, it drops out of the 
problem. The only relevant scale is the density and as a result the gas has the same scaling behavior as an ideal 
Fermi gas. This allows one to fix immediately the leading power 71 = 2/3, although the next term in the expansion 
is not analytically known. It is commonly proposed to be linear in the parameter x = — 1 / fcpa s [T5L liH IrH Il8|. 
Thus the expansion of the chemical potential can be written as fi(x) = (£ + Sx)h 2 kp/2m^ where the dependence on 
the density is accounted for by the Fermi momentum kp — (S^n) 1 ^ 3 . We find that the predictions of the sum-rule 
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approach (see Table ITT)l coincide in the unitary regime with the results of ^||, where the authors start directly from 
the hydrodynamic equations and solve them for this particular choice of the equation of state. The second row of the 
Table |n] shows the predictions in the Bardeen-Cooper-Schrieffer (BCS) limit a s — > —0 (A^d — * 0), where fermions 
experience a weak attraction and the equation of state can be found perturbatively starting from a weakly interacting 
Fermi gas. The leading term corresponds to an ideal Fermi gas and scales as n 2 / 3 with the density. This leads to the 
same leading term in the frequencies of the collective oscillations as at unitarity, although the next term [Ifi, l20j is 
different and arises from the interparticle interaction of a normal Fermi gas. 

The situation in the Bose-Einstein condensation (BEC) regime a s — > +0 (Ajd —> 0) is extremely interesting. In 
this limit a two-body bound state (molecules) of atoms of different spin is formed. The main contribution to the 
energy comes from the negative binding energy of molecules. As the binding energy depends only on the interaction 
potential and is independent of the density limit, it docs not contribute to the collective excitations and this term 
is omitted in Table [H] The first leading term depending on the density corresponds to a dilute non-ideal Bose gas 
of molecules, which play the role of composite bosons with mass 2m and density n/2. The s-wave scattering length 
a m of such composite bosons was found directly by solving the four-body scattering problem [2l| and was estimated 
to be a m — 0.6 a s . It should be noted that this result differs from the commonly used mean-field BCS theory which 
gives a m = 2 a s (see, for example, |22l l23l l24j1. The nature of the the next term in the equation of state is still an 
open question and only a precise experimental measurement can give a definitive answer. This term might have a 
bosonic nature and correspond to the beyond mean-field correction |25| to the equation of state of composite bosons 
as proposed in [2(| (see, also, H3)- Alternatively, the next correction might be of a fermionic nature as argued in |28j . 
An ab initio quantum Monte Carlo study of the many-body problem made it possible to obtain the equation of state 
numerically [2S|. Firstly, it was shown that the leading term has a mean- field nature and agrees with a molecule- 
molecule scattering length equal to a m — 0.6 a s (very recently this result was also recovered from a diagrammatic 
approach J3(J]). Secondly, the Monte Carlo study suggests the presence of the bosonic beyond mean-field correction. 
As inferred from Table ITT1 and shown in Fig.^this would mean that the frequencies of the collective oscillation should 
increase above the mean-field value by making the interaction between molecules stronger. To date there has been 
no experimental corroboration of this result. In the available experimental work [T^ | the frequency decreases, but the 
deep BEC regime has not yet been achieved. A more precise experimental study should be performed, in which case 
Eq. I|16l) can be used to fit the frequencies of the collective oscillations and thus provide insight into the nature of the 
equation of state. 

We provide a relationship between the characteristic LDA parameter A3D and the value of the Fermi momentum in 
the center of the trap kF,max, which is a natural quantity for a fermionic system. The leading term (which is sufficient 
for the applicability of Eq. I|16|l ) is given by 
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TABLE II: Summary for three-dimensional systems. The first column labels the physical system under considerations: limits 
for the two-component Fermi gas in the BEC-BCS crossover and weakly-interacting system of bosons. Columns 2-5: coefficients 
of the expansion of the equation of state defined as in Eq. ©. Columns 6-7: frequencies of radial and axial oscillations in 
a "cigar" -shaped trap (note that in a 3D system the parameter lib! is 3 = Q^/uj'^). Column 8: frequency of the breathing 
mode in a spherical trap. The notation a m in the BEC regime denotes the molecule-molecule scattering length in units of the 
atomic scattering length. In the same regime we omit the contribution to the chemical potential from the binding energy, as it 
is independent of the density and, thus, is not important for the LDA. 



Low-dimensional geometries are advantageous for the investigation of the energetic properties of a gas as the 
reduced phase space enhances the role of interactions. As a result, the frequencies of collective oscillations vary over 



7 



2.54 



2.52 



2.50 
2.48 
2.46 
2.44 
2.42 



2.40 

0.0 0.4 0.8 1.2 1.6 

x . -1/2 

atan A 3D 



FIG. 1: Square of the axial mode frequency in "cigar" -shaped trap for the Monte Carlo (thick solid line) |29| and the mean-field 
BCS (thin line) [23, E3, equation of states as a function of A313. Experimental results (symbols) are taken from |l3j. The 
dashed lines show the expansions at unitarity (short-dashed line) and in the BEC regime (long-dashed line) taken from Table HT1 

a larger range, thus facilitating possible experimental measurements. We consider a quasi-two-dimensional system 
of bosons in the regime of small densities. We assume that the gas is weakly bound in two-directions (x,y) and 
is subjected to a tight confinement in z direction. The density profile in a uniform system can be approximated 
as n(x,y,z) — n2D exp(— z 2 ) / ^Tra 2 sc , where n2D is the two-dimensional density. Locally the equation of a rarefied 
three-dimensional bose gas can be applied (see Table [HJ and the energy of the gas is estimated by integrating out 
the variable corresponding to the tight confinement ^213(^213) = / E^jj(n(x,y, z))n(x,y, z)dz/ J n{x,y, z)dz. The 
chemical potential, obtained by differentiating the energy, is reported in Table II I II We calculate frequencies of 
collective oscillations by exploiting Eqs. H16ll7fl . The results for a spherical and very elongated traps are presented 
in Table ITTT1 We note that the quasi two-dimensional equation of state is valid when the scattering length a s is small 
compared to the oscillator length of the tight confinement, so that the scattering process is still three-dimensional. 
In a "pure" two-dimensional system the equation of state contains logarithmic terms |3l| and our approach does not 
apply. 
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Ci 


71 


C 2 


72 




0.1/ Ulz 




Q2D Bose gas 


2v / 2lr" 


1 


64 


1/2 


2 5 / 8 I92A^ 2 


8 2 5/8 256A^ 2 


2 5 / 8 384A^ 2 
4 H = — — 




35\/57r 9 / 8 


3 105x/57r 9 / 8 


35\/57r 9 / 8 



TABLE III: Equation of state of a quasi-two-dimensional bose gas in the weakly interacting regime. The expansion parameter 
is ii2Dffl 3 / fl o»c- ^ i s assumed that a <C a osc . The meaning of the columns is the same as in Table UTI 

In Table llVl we summarize the examples of bosonic as well as fermionic systems in (quasi-) one-dimensional geometry. 
The Lieb-Liniger model |32| describes a gas of repulsive bosons interacting with a contact (5-potential. The equation 
of state of such a gas can be expanded in the weakly interacting limit Aid = ^ a iD/ a ho ~ * 00 ( one notes that in a 
waveguide small values of the three-dimensional s-wave scattering length a^D correspond to small values of the one- 
dimensional s-wave scattering length a\u as in this limit auj = —a,j_/a z ). The leading term in the chemical potential 
jj, is linear in the density and corresponds to the mean-field Gross-Pitaevskii theory. Zero point oscillations of quasi- 
particles contribute to the next term and in a three-dimensional system corresponds to the Lee- Yang correction |2C| . 
We applied the Bogoliubov theory to a one-dimensional system and found that the corresponding term coincides 
with the second term of the expansion of Lie-Liniger equation of state. In the opposite limit of strong correlations 
Aid — * (Tonks-Girardeau regime) the leading term in /1 equals the Fermi energy of ideal spinless fermions due to a 
very peculiar property of a Tonks-Girardeau gas experiencing "fermionization" , as the wavefunction of impenetrable 
point-like bosons can be mapped directly onto a wavefunction of a fermionic system |3^ | . Another one-dimensional 
system, the Calogero-Sutherland gas j3J, |35| has a functionally similar equation of state (see Table IIVI parameter A 
is related to the strength of the l/|.z| 3 interaction). In order to understand the physical meaning of the subleading 
correction we will consider for a moment a system of hard-rods, i.e. impenetrable bosons of a size |aiu|. Its equation 
of state can be obtained from the energy of a Tonks-Girardeau gas by reducing the density by the excluded volume 
n — » n — N\clid\- The leading term is again given by the energy of an ideal Fermi gas of spinless particles and the next 
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correction is given in Table llVl It was shown in [3g that a lowest gas- like state of attractive bosons (am > 0) in the 
limit of small (super- Tonks gas) has the same leading term in the expansion of the equation of state as a system 
of hard-rods. The reason is that the node of a two-body low-energy scattering solution lies at a\o and for larger 
distances is essentially the same as in a system of hard-rods. Instead, for a repulsive gas with a contact interaction in 
the vicinity of the Tonks-Girardeau limit, the analytic continuation of the scattering solution has a node at — |<xix> ] - 
One finds that the subleading term in the equation of state is equal to the "excluded volume" correction, but differs 
in sign (see Table ITvTl . Of course, the discussed terms in the chemical potential of the Lieb-Liniger gas coincides with 
the ones found by solving iteratively the Bethe ansatz integral equations. Still we find it important to understand the 
physical origin of the terms of the expansion. 
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Lieb-Liniger: strong interaction 
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Attractive Fermi gas: strong interaction 
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1 
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Attractive Fermi gas: weak interaction 


7T 2 /8 


2 


-8/tt 2 


-1 


4 + M 1/2 


Repulsive Fermi gas: strong interaction 


7T 2 /2 


2 


-81n(2)/3 


1 


A 128v / 21n2 Al/2 


Repulsive Fermi gas: weak interaction 


tv 2 /8 


2 


8/^ 2 


-1 


4_^L/A 1/2 


Gas of Hard-Rods 


7T 2 /2 


2 


8/3 


1 


128V2 1/2 

4+ 1.^2 ID 


Calogero- Sutherland 


AV/2 


2 


8/3 


1 
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TABLE IV: Summary of some of the one-dimensional models where the expansion of the equation of state is known. The first 
column names the considered systems. The coefficients of the expansion are given in columns 2-5. The last column gives the 
predictions 1181 for the oscillation frequencies. The parameter Aid is defined by 0. Note that the presence of a term in the 
chemical potential independent of the density (e.g. binding energy of a molecule) does not modify the frequencies of oscillations 
and is ignored. 

Considering a system of one-dimensional fermions and limiting ourselves to the s-wave scattering we immediately 
find that atoms of the same spin do not interact due to the Pauli exclusion principle. Instead, fermions of different 
spin are allowed to interact and the interaction can be either attractive or re puls ive. The equation of state can be 
found by solving the integral Bethe ansatz equations as discussed for attractive |37ll38| and repulsive [3^ interactions. 
The numerical solution for the LDA was carried out in |4Cj and is compared here to the series expansions obtained 
from Eq. 1)18(1. Again, there is interesting physics behind the expansion J7J. In the weakly interacting limit (both 
attraction and repulsion) the leading term equals the energy of an ideal two-component Fermi gas. The next correction 
is linear in density and describe the mean-field interactions between particles of different spin. In the limit of strong 
repulsion the main contribution to fi coincides with the chemical potential of a one-component ideal ID Fermi gas 
with N atoms, as the strong repulsion between atoms of different spin plays the role of an effective Pauli principle 
|4l| . An attractive interaction in one-dimension immediately leads to the formation of a bound state (molecule) . As 
the binding energy is independent of the density and does not contribute to the collective oscillations, this term is 
omitted in Table ITT1 In the low-density limit the internal structure of a molecule of two fermions can be neglected and 
be considered as a composite boson. Since in this limit the interaction between such composite bosons is repulsive 
and strong (Tonks-Girardeau limit of N/2 bosons of mass 2m), the system can be mapped to N/2 spinless fermions. 
So the leading behavior is again of a fermionic nature. The next term was obtained correctly for the first time in 
[Hill and is reported in Table llVl 

In Figs. |21 El we visualize how the terms of the EOS expansion show up in the frequencies of collective modes. The 
Fig. |21 shows the frequencies for bosonic systems, although the Fig. |3 presents results for systems of two-component 
fermions. We compare the predictions of the expansion, Eq. (|18fl . (dashed lines) with the numerical solution for the 
exact equation of state (solid lines). As it can be seen from the Figures, the region where the perturbative results 
hold is quite large, making possible the extraction of the expansion terms of the equation of state from the frequencies 
of collective oscillations. The accuracy is better for the expansions at low density, since if the expansion holds for 
the center of the trap it will certainly be valid close to the borders of the trap where the density is smaller (see, e.g. 
Fig.EJ. 
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FIG. 2: Frequency of the breathing mode in a number of one-dimensional bosonic systems. The frequencies for the Lieb-Liniger 
model were obtained numerically in ;9]. The collective excitations of a gas of hard-rods were discussed in [3(J in their relation 
to the excitations of a gas in the super- Tonks regime. The dashed lines correspond to the expressions given in Table HVI 
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FIG. 3: Frequencies of the breathing mode of a two-component Fermi gas in a quasi one-dimensional geometry. The solid line 
shows results obtained by solving numerically the LDA with an exact equation of state [4(1 ]. The dashed lines show the small- 
and large- density expansions from Table HV1 The upper curve corresponds to attraction between atoms of opposite spin, while 
the lower curve corresponds to repulsion. 

VI. CONCLUSIONS 

It quite often happens that a perturbation theory for a homogeneous system can be developed analytically, although 
no solution for a system in an external confinement is known. If the number of particles in a trapped system is large 
enough (and this requirement often can be easily fulfilled in experiments) the local density approximation can be 
applied. We assume that the series expansion of the dependence of the chemical potential on density is known for a 
homogeneous system and we find the properties of a gas confined in an anisotropic trap. We consider three-, two-, 
and one- dimensional systems. We obtain the energetic properties such as the chemical potential, total energy of the 
gas, potential energy of the confinement. The calculated release energy can be accessed in expansion experiments. 
Spatial properties such as the density profile and size of the cloud are predicted and can be measured by imaging 
techniques in current experiments. 

The frequencies of the collective oscillations in a harmonic trap are predicted from a hydrodynamic theory. We 
apply a scaling approach to obtain frequencies in three-, two-, and one- dimensional systems in traps with arbitrary 
aspect ratio luj_/lu z . In spherical, "cigar"- and "pancake"- shaped geometries the expressions are simplified. In a 
spherical trap an equivalent method, the sum-rule approach, can be conveniently applied and leads to analogous 
expressions. Exploiting the obtained expressions and using our solution to the LDA problem we relate the frequencies 
of the breathing modes to the expansion of the equation of state of a homogeneous systems in an explicit way. 

We argue that the obtained results are highly relevant for the experimental study of the equations of state of 
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trapped gases. The expansion of the equation of state can be found by a best fit of our formulas to experimentally 
measured quantities (for example, collective oscillation frequencies). This would provide an important insight into 
the physics of the systems under the study. As an example, there is an open and intriguing question of what will be 
the nature of the beyond-mean field term in the equation of state of the recently observed condensation of molecules. 
Will it have a bosonic or fermionic nature? Both proposals can be found in the literature and only a high precision 
experiment can provide a definitive result. 

We apply the derived formalism to a number of three-dimensional systems (two-component Fermi gas in the BEC- 
BCS crossover, dilute Bose gas), quasi two-dimensional bose gas, quasi one-dimensional systems (Lieb-Liniger gas, 
gas of hard-rods, two-component fermionic gas, Calogero-Suthcrlanf gas). A comparison with numerical solutions 
is provided. Finally we provide an physical interpretation of the expansion terms in a number of one-dimensional 
systems. 

This work was supported by Ministero dellTstruzionc, deU'Universita e della Ricerca (MIUR). Useful discussions 
with S. Giorgini and S. Stringari are gratefully acknowledged. We thank B. Jackson for reading the manuscript. 
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